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Nodal lines (phase singularities, optical vortices) are the generic interference 
fringes of complex scalar waves. Here, an exact complex solution of the time in- 
dependent wave equation (Helmholtz equation) is considered, possessing nodal lines 
which are braided in the form of a borromean, or pig-tail braid. The braid field is a 
superposition of counterpropagating, counterrotating, non-coaxial third order Bessel 
beams, and a plane wave whose propagation is perpendicular to that of the beams. 
The construction is structurally stable, and can be generalized to a limited class of 
other braids. 

PACS numbers: 42.25.Hz, 03.65.Vf, 02.10.Kn 

In three dimensional fields of interfering complex scalar waves, perfect destructive in- 
terference, that is, where the total wave amplitude is zero, generically occurs along lines: 
nodal lines. They are also called phase singularities, since the phase is undefined when 
the amplitude is zero, or wave dislocations 0, Q]. Nodal lines in waves are much studied, 
particularly in solutions of the time independent (Helmholtz) wave equation, in which they 
are stationary features of the interference pattern. They appear in a variety of physical con- 
texts, including quantized vortex lines in solutions of the Schrodinger equation in quantum 
mechanics Isll, and optical vortices in scalar components of monochromatic electromagnetic 
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fields 12, li, where t heir presence in optical beams often indicates the presence of orbital 
angular momentum ^. Generally, in interfering spatial fields, there are very many nodal 
lines, tangled in a complicated way [^, reminiscent of the tangling of quantized vortices in 
media such as turbulent superfluid helium 

Here, I shall demonstrate that a particular complex wave superposition (involving two 
Bessel beams, defined in equation (0) below, and a plane wave) contains braided nodal 
lines. As shown in figure ^ the field has three phase singularities whose trajectory, as z 
increases, traces out a figure-8 in the transverse xy-plane, forming the well-known borromean 
braid (also called the pigtail braid, or plait). This crafting of waves with a braided nodal 



FIG. 1: Braided nodal lines in the field described by equations ©, ©• The black curves represent 
the X, y trajectories of the braid strands, and the structure is plotted to scale. Three z-periods are 
plotted, between z = and Qir/kz- 



configuration complements wave fields containing knotted and linked nodal lines in free space 
0,0] and in the quantum mechanical electronic states of hydrogen [l^. Braiding of nodal 



lines in light also complements Ref. in which two laser beams in a plasma were twisted 
into a double helix; the topology here is more complicated (the transverse trajectory being 
a figure-8, rather than a circle), and occurs in linear wave superpositions in free space. The 
topology of braided nodal lines is reminiscent of braided vortices in topological fmagneto-) 
hydrodynamics 12], and braided vortex trajectories in hamiltonian mechanics 13 1. 

It is easy mathematically to find complex scalar fields which contain various arrangements 
of nodal lines Q], as products of separate complex functions for each nodal line. For each 
such line, two independent real functions ^(r),?7(r), are found whose zero contours cross on 
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the desired line; the function ^ + iri has a nodal line in this position. For example, the field 
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V'Liss = - + 2^j73) + Ky - sin 2(2; + 2ttj/3)) (1) 
3=0 

has three nodal lines, and as z increases, each executes a simple figure-8 Lissajous trajectory 
in the xy-plane. The total field is the product of three complex fields; it has the zeros of 
each of its separate product terms, and the nodal structure of the total field is braided. 

However, such fields are not, in general, physically realisable as waves. That is, the field 
dH) does not satisfy the time independent Helmholtz equation 

VV + = 0. (2) 

Nodal lines are usually studied in solutions of this equation because they are stationary, 
even if the field is changing periodically in time (as in monochromatic optical fields). In 
the following, the wavenumber k will be set to 1. The most familiar set of solutions to (j21) 
are plane waves exp(ik ■ r) where |k| = 1; these do not themselves have phase singularities, 
although superpositions of them can. 

Associated with phase singularities is a topological number, the topological current (or 
strength). This is an integer, equal to the total phase change around the singularity divided 
by 2tt. The right-hand sense of phase increase around the line endows it with a direction; if 
the nodal line crosses a plane at a nodal point, the associated topological charge is positive 



if the sense of the line is out of the plane, negative if into 



1^]. A convenient set of solutions 



of O), which contain a strength m phase singularity up the 2;-axis are the Bessel beams 



(, . \ m/2 
^ _ . ] exp(iA;^2;). (3) 



The factor ^Jl — k1 in the argument of the order m Bessel function is the transverse 
component of the wavevector (fixed by fc^, since the total wavenumber is 1); the Bessel 
function has a nodal line along x = ?/ = 0. The ratio oi x + iy and x — iy io the m/2 power 
in (j2)) gives an order m phase singularity: if m > 0, the direction of the topological current 
is in the -|-z-direction, if m < 0, it is in the — 2;-direction. The beam itself propagates in ±0 
according to the sign of kz. An advantage of using these nondiffracting singular beams here 
(instead of, for example, Laguerre- Gauss beams) is that their structure in the 2;-direction is 
periodic with period ^Txjkz- 
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High order singularities, when |m| > 1, are not structurally stable. If a Bessel beam Q is 
perturbed by a function which breaks its azimuthal phase pattern, the singularity generically 
unfolds to |m| unit strength nodes, wound as a helix; the singularities are convected by 
the twisted phase structure of the beam. Unlike hydrodynamic vortices, phase singularity 
vortices do not repel, but may cross, in a manner equivalent to reconnection of vortices in 
superfluids The crossings themselves are not stable, and, on perturbation, break in 

either of the two ways that conserve topological current j9|. 

The knot construction of Ref. 0| employed a superposition of n{n + l)/2 parallel, coaxial 
Bessel beams of the same order m, with different kz and coefficients, in order to make a 
strength n nodal loop at fixed radius from the beam axis, threaded by the strength m 
axial nodal line. On perturbation by a Jq beam, the loop unfolds to give an (m, n) torus 
knot or link, which winds around the z-axis m times as it winds around the original loop 
position n times. The braid construction here employs a similar technique, by superposing 
parallel Bessel beams, in this case counterpropagating and non-coaxial, to get a degenerate 
configuration, which unfolds to a braid after perturbation by an appropriately chosen plane 
wave. 

The solution ipbrmd of (0), containing the braided nodal lines in figure [TJ is given by 

'4'hri>.id{x,y,z) = ip3{kz,x + d,y,z) + 'ip3{-kz,x- d,y,z) + ieexp{-iy); (4) 

where the parameters e, d and kz have values 

e = 0.05; = V37r/2 = 2.721 ... ; A;^ = (1 - jVStt^)^/^ = 0.772 ... , (5) 

and j = 4.201 . . . , the first positive zero of J3. The reason for the particular choices (0) is 
explained later. The first two terms of Q are order +3 Bessel beams of the form Q, the 
first propagating in +z and centered on the line x = —d,y = 0; the second propagates in 
—z, and is centered on x = +d, y = 0; the order 3 topological current in both beams is in the 
+z direction. The third summand is a plane wave propagating in the — y-direction. Since 
these three waves propagate in very different directions {+z, —z and —y), the braid function 
(0} is far from the paraxial regime. 

I begin by describing the nodal configuration of V'braid without the perturbing plane wave 
ieexp(— iy), that is, 

An{x, y, z) = ^jJsikz, x + d,y,z) + ^psi-kz, x-d,y,z). (6) 
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FIG. 2: Configuration of tlie unperturbed braid field V'un (© (that is, the braid field ^braid @ 
with the plane wave amplitude e = 0). The two horizontal planes are at z = and 27: /kz (one 
z-period). The viewpoint is from the — x-direction. 
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FIG. 3: Phase pattern of the unperturbed configuration Q in the plane z = 0. There are phase 
singularities at (x,y) = (0,0) and (0, ib37r/2). 

The configuration of nodal lines of © in the vicinity of the 2;-axis is shown in figure El It 
shows a series of tilted singularity rings encircling lines confined to the a; = plane. The 
rings meet the threading lines at nongeneric connections, at 

a; = 0, y = ±V?>d = ±37r/2; z = 0, n/k^ mod 27T/kz. (7) 
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The phase structure in the z = plane is shown in figure El There is one singularity at the 
origin (this is one of the threading curves in figure Ej), and two others at x = 0,?/ = ±37r/2. 
These latter nodes are nongeneric, since not only are they zeros of the field, but also places 
where the x-derivative vanishes. This refiects the fact that in three dimensions, they are 
crossings ((Zj), shown in figure El The position of phases around the nongeneric singularities 
is the same for each (e.g. =F7r/2 in ±y directions, phase increases in the same sense around 
each) . 

The configuration of singularities of figure |21 appears characteristic of the nodal pattern 
of superposed non-coaxial Bessel beams which are counterrotating (i.e. counterpropagating 
with parallel topological currents), and contains the basic structure of the three strand 
figure-8 xy-trajectory desired for the braid. However, it is only with the particular choices 
of (0) that the crossings occur at precisely the same z value, which is convenient for the 
present purpose; increasing d slightly makes the two crossings migrate away from the z = 
plane (and each other), d and kz were chosen in so the crossings occur in the same planes 
z = 0,Tr/kz (mod 2n/kz), and with a spacing of Stt in the y-coordinate. The crossings ((Zj) 
are broken when the perturbing plane wave ieexp{—iy) is included in 

The pair of crossings at 2; = in ((7|) are separated by a distance of Stt; that is, one and a 
half wavelengths of the plane wave. Thus, the phases of the plane wave at the positions of 
the two z = crossings of tp^^ differ by tt. This causes the two crossings to break in opposite 
senses (figure E}, and the two zeros in the z = plane of the braid field are perturbed in 
opposite directions, as in figure El Moreover, at z = ir/kz (half the z-period), all of the 
phases of ip^n are reversed. The plane wave, however, is not changed since its propagation 
is in the ^/-direction alone; this implies the sense in which the crossing is broken by the 
plane wave is opposite to that of the corresponding crossing at 2; = 0; the sense of crossing 
breakings alternates as z increases. The full set of nodal lines near the 2-axis, shown in figure 
El includes the braid, shown in black; there are other nodal lines, shown in grey, which have 
broken away from the configuration of figure El and do not take part in the braiding. 

The braid described here is structurally stable; if the numerical values of the parameters 
in (0) are changed slightly, the braid remains topologically unchanged, although past certain 
critical values of kj,, d, and e, it will be dissolved through various reconnection events, as 

n 

with nodal knots 9]. Structural stability means that parameters such as d and kz need not 
be precisely controlled in an experimental or numerical realisation. 
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FIG. 4: Perturbed braid configuration, where the field of figure has been perturbed by the plane 
wave 0.05iexp(— iy). The black curves form the braid plotted in figure^ the grey curves represent 
other nodal lines. As with figure [21 the distance between the two horizontal planes is one z-period 




X 



FIG. 5: Phase pattern of the full braid field V'braid in the plane z = 0. The figure-8 x,y trajectory 
of the braiding singularities is included. The singularities depicted as grey in figure 0] do not cross 
the z = plane, and so are not seen here. 
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The key feature of the braid reahsation is that the trajectory of the strands of a borromean 
braid is topologically a figure-8, which has been achieved by wave dislocations in superposed 
off-axis counterrotating beams. The parameters d and kz were chosen in ^ so that the 
perturbing plane wave would have opposite phases at the two reconnection points with 
z = 0; the reconnections were separated by Svr. However, this argument works for any 
separation 7r(n + 1/2), for n = 1, 2, ... , for which the values of d and kz in (0) should be 
replaced by 

d = n{n + 1/2) /VS; fc, = (1 - 3f/n^i2n + 1)^'. (8) 

In (jSj), n = 1 (if n = 0, then k"^ < 0, and the wave is not periodic in z). In this case, e can 
be arbitrary but small; the braid structure is more sensitive to the size of perturbation for 
higher n. V^braid can straightforwardly be generalised to give other m-stranded braids whose 
transverse trajectory is a figure-8, by using Bessel beams of odd order m > 3, as follows. In 
the unperturbed situation, with 3 replaced by m, there is a series of m zero points in the 
z = plane on the x = line; kz and d are chosen such that the two furthest from the origin 
are ai y = ±TT{n + 1/2), and that the x-derivative vanishes there also. If m is even, the 
extremal y zeros must be separated by an integer number of wavelengths, and the perturbed 
pattern, instead of giving a braid, gives a pair of m/2-stranded helices. The construction 
cannot be generalized in any obvious way to give other, more general braids. 

Bessel beams and plane waves are often used in optics, but it is not entirely clear what 
sort of quantum wave would correspond to (jH). Would it therefore be possible to create a 
close enough approximation to V'braid in an optical experiment? This leads to the following 
observations. Electromagnetic waves are of course vector fields, and, in general, phase 
singularities do not occur in vector waves, especially in the nonparaxial regime [2^; nodes 
generically only occur in scalar components of the wave. In the present case, if all three wave 
beams in are linearly polarized in the x-direction, which is possible since the separate 
wave beams are propagating in the ±z- and — y-directions only, the braided nodal lines 
appear in the x-component of the total field (there will also be a longitudinal 2;-component 
arising from the Bessel beams). A more serious problem in optical realization is the departure 
from paraxiality; not only in the interference of coherent counterpropagating beams, but in 
the fact that in the Bessel beams kz ~ 0.8, meaning that the divergence of the beam is large. 
Moreover, the braid period is on the scale of the beam period itself, which is extremely small 



at optical frequencies - unlike the knotted singularity construction of Ref. 



, the absolute 
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scale of the braid structure is fixed. However, it is possible that use of a higher n as in 
(jHI) may lead to a more experimentally realistic beam. It may be possible to use fields in a 
different part of the electromagnetic spectrum, such as radio waves, for which Bessel beams 
can be synthesized and studied 21 1. 

The present work shows how potentially realizable beams may be superposed to form 
a field containing structurally stable, braided nodal lines. However, this realization of a 
braid is not necessarily the only possible one in interfering waves. For instance, it would 
be interesting to investigate whether periodic braiding structure is possible in the space 
above a two dimensional grating. Such a structure would have the desired 2;-periodicity, 
although it is not composed of counterpropagating waves as used here. Braids, like knots, 
are particular cases of a wide range of topologically interesting phenomena that may occur 
in three dimensional wave fields, and one may ask whether other nodal configurations, such 
as chains, are possible in beam superpositions. The existence of braided nodal lines widens 
the topological possibilities of singular optics and quantum vortices, and provides a technical 
challenge for experiment. 

I am grateful to John Hannay and Michael Berry for discussions. The figures were 
generated using Mathematica; figure^ with Mitchell Berger's tuba package. This work was 
supported by the Lever hulme Trust. 
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